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ABSTRACT 

We have analyzed the kinematical parameters of CirX-1 to constrain the nature of its 
companion star, the eccentricity of the binary and the pre-supernova parameter space. 
We argue that the companion is most likely to be a low-mass 2.0 M©) unevolved 
star and that the eccentricity of the orbit is 0.94 ± 0.04. We have evaluated the dy- 
namical effects of the supernova explosion and we find it must have been asymmetric. 
On average, we find that a kick of '-^^740 km s^^ is needed to account for the recently 
measured radial velocity of -1-430 km s^^ (Johnston, Fender & Wu) for this extreme 
system. The corresponding minimum kick velocity is ~500 km s~^. This is the largest 
kick needed to explain the motion of any observed binary system. If CirX-1 is asso- 
ciated with the supernova remnant G321.9-0.3 then we find a limiting minimum age 
of this remnant of ^ 60 000 yr. Furthermore, we predict that the companion star has 
lost ~ 10 % of its mass as a result of stripping and ablation from the impact of the 
supernova shell shortly after the explosion. 
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1 INTRODUCTION 

CirX-1 is a unique binary in many aspects. Despite inten- 
sive X-ray observations since its discovery in the early 1970s 
very little is known about its stellar components besides an 
orbital period of 16.6 days (Kaluzienski et al. 1976). It is 
still uncertain whether the binary is a high-mass (HMXB) 
or low-mass X-ray binary (LMXB) . The discovery of Type I 
X-ray bursts (Tennant et al. 1986) demonstrates that the ac- 
creting compact object must be a neutron star. In a recent 
paper Johnston, Fender & Wu (1999) presented new optical 
and infrared observations of CirX-1 which show asymmetric 
emission lines. Combined with twenty years of archival Hq 
line emission data, they interpreted that the narrow com- 
ponents of the lines imply a radial velocity of -1-430 km s~^ 
for the system (including a small correction for the Galactic 
rotation). This is the highest velocity known for any X-ray 
binary. 

CirX-1 is located 25' from the centre of the supernova rem- 
nant G321.9-0.3 and is apparently connected by a radio neb- 
ula (Haynes et al. 1986). Furthermore, there is some evi- 
dence for an association from observations of arcmin-scale 
coUimated structures within the nebula (Stewart et al. 1993) 
which are likely to originate from an arcsec-scale jet which 
has been observed (Fender et al. 1998) to be aligned with 
these larger structures. If the suggested association is cor- 
rect then CirX-1 is a young (< 10"" yr) runaway system 
and, in combination with the recent estimated distance to 



the remnant of 5.5 kpc (Case & Bhattacharya 1998), the 
inferred minimum transverse velocity (390 km s"'^) yields a 
resulting 3D space velocity of >580 km s~^. In this Letter 
we investigate what can be learned about CirX-1 and the 
effects of the supernova explosion which gave rise to this 
high runaway velocity and orbital period of 16.6 days. 



2 DYNAMICAL EFFECTS OF AN 

ASYMMETRIC SUPERNOVA EXPLOSION 

We use the analytical formulae by Tauris & Takens (1998) 
to calculate the dynamical effects of an asymmetric super- 
nova (SN) in CirX-1. These formulae also include the effect 
of shell impact on the companion star. We assume that the 
exploding He-star with mass Mhe (the progenitor of the neu- 
tron star with mass Mns) is at the origin of the cartesian 
coordinate system which we shall use in our description. 
The positive z-axis points in the direction of the pre-SN or- 
bital angular momentum. The positive y-axis points towards 
the companion star (with mass M2) and the positive x-axis 
points in the direction of v, which is the pre-SN relative ve- 
locity vector of the He-star with respect to the companion 
star. We assume the pre-SN orbit is circular and ro is the 
pre-SN separation between the two stellar components. The 
systemic velocity of a binary which survives the SN is: 



VApf+Apf+Api/(MNS + M2f) 



(1) 
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where the change in momentum resulting from the SN is; 



MnsM2 - MHeM2f 



Mho + M2 

Apy = A/ns w sin cos ip + M2{ vi 
Apz = A/ns w sin i9 sin ip 



V + Mns w cos 1? 



(2) 



A kick with magnitude w is imparted to the neutron star 
(NS) by an asymmetric explosion. The pre-SN relative or- 
bital speed IS V = G{Mnc + M2)/ro and the angle be- 
tween V and w is i9. The second position angle ip of w, is 
taken such that Wy = w sin i? cos ip. Finally, Vim is an effective 
speed taking into account the combined effects of incident 
shell momentum and the subsequent momentum resulting 
from mass loss due to stripping and ablation of stellar ma- 
terial from the surface layers heated by the passing shock 
wave, ft is given by (Wheeler, Lecar & McKee 1975; see also 
Fryxell & Arnett 1981): 



Vim = rj Uojcct 



/ -R2a:crit \ / A-/shcll \ 

V 2ro J \ M2 J 



1 + ln(2 

I' eject / ''^ei 



(3) 



(-Pstrip ~t- -/^ablate) 

where R2 is the initial radius of the companion star, 
'^ojcct is the speed of the material ejected in the SN, 
Vesc = \/2GM2{xcvit) / R2 Xc-At is the escape velocity from 
the stripped companion star and the parameter ^^-crit IS a 
critical fraction of the radius outside of which the total 
mass fraction, Fstrip is stripped and inside of which a cer- 
tain fraction, Fabiatc is ablated. Thus after the shell im- 
pact the new mass of the companion is given by: M2i ~ 
M2 (1 — (-Fstrip + ^kbiate))- A correctiou to the assumption of 
plane parallel layers is expressed by -q. We assume rj — 0.5. 
Wheeler, Lecar & McKee (1975) have defined a parameter, 
* in order to evaluate Xcm and F = Fstrip + -Fabiate (their 
tabulation is given in Table 1): 



* = (l^) (1f)(r-') 

In order to calculate the possible values of 'I' for the case 
of CirX-1 we assume that a main-sequence companion star 
can be approximated by a polytropic model with index n—3. 
In this case one has: M2{xcTit) / R2 Xait — 5/3 MqRq^ for 
0.2 < Xcrit < 0.7 (to an accuracy within 15 %). Furthermore 
we assume Sojcct = 1/2 Afaheii «cjcct = 1-0 x 10^^ erg. 
The absolute minimum ZAMS mass for producing a NS is 
~ 8 Mq . Such an early-type star evolves on a short timescale, 
and hence at the time of the SN the maximum age of the 
companion is ~40 Myr. The companion star was less mas- 
sive, and is therefore most likely to be unevolved. We fitted 
a mass-radius relation for ZAMS stars using the stellar evo- 
lution models of Pols et al. (1998): 



R2/R(i 



0.86 (M2/M0 
1.15 (Ma/Af© 



|1.21 
vO.50 



M2 < 1.5 Mq 

M2 > 1.5 Afo 



(5) 



assuming X=0.70, Z=0.02 and a convective mixing-length 
parameter a = 2.0. 

The post-SN separation and eccentricity are given by: 
ro 



" 2-e e = vi + (?-2)(g + i) 

where (see Tauris & Takens 1998): 

+ + + 2w {v cos ■& — Him siu cos ip) 



(6) 



(7) 



Table 1. Parameters as a function of ^ (see text) for a companion 
with the structure of a polytrope with index n=3. 
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^crit 


-^strip 


-^ablate 
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0.29 
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0.26 


0.43 


0.34 


0.77 
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0.29 


0.35 


0.34 


0.69 
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0.37 


0.19 


0.28 


0.47 
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0.44 


0.10 


0.19 


0.29 


0.4 


0.52 


0.04 


0.11 


0.15 


0.2 


0.58 


0.02 


0.07 


0.09 


0.1 


0.64 


0.01 


0.04 


0.05 


0.01 


0.78 


0.001 


0.006 


0.01 


0.001 


0.87 


~ 


~ 


~ 



This table was taken from Wheeler, Lecar & McKee (1975). 



{w sin 1} cos p - 



(8) 



and m = {Mns + M2i)/{Muc + M2). 

Prior to the SN we require that the companion star is able to 
fit inside its Roche-lobe, and after the SN it is required that 
the companion star must fill its Roche-lobe near the post-SN 
periastron, a (1 — e) in order to explain the observed X-ray 
emission (which is caused by mass accreted onto the neutron 
star from a disk which has to be fed by stellar material). 
Combining our mass-radius relation with Eggleton's (1983) 
expression for the Roche-lobe radius, Rl yields the following 
criterion for the radius of the companion star: 

RUqo)Uro >R2> R2f > i?L(gf)l,,=,(i_,) (9) 

where Ri,{q) = r 0.49 g^/^ (0.6 g^^^^ + ln(l + q^^^))~^ and 
the mass ratios before and after the SN are given by: 
go = M2/M-RC and qf = Af2f/AfNS, respectively. 



3 SIMULATION RESULTS AND DISCUSSION 

We have performed Monte Carlo simulations of the SN ef- 
fects on a large number of binaries to determine the mass 
of the companion star, the eccentricity of the present or- 
bit and the pre-SN orbital parameter space. The constraints 
which must be fulfilled by any post-SN binary in order to be 
consistent with present observations of CirX-1, are: I) the 
companion star must fill its Roche-lobe at periastron, II) 
the system must have a high velocity of at least 430 km s~^ 
(Johnston, Fender & Wu 1999) and III) the system must 
have an orbital period between 14 and 18 days. For CirX-1 
Poih ~ 16.6 days, but tidal effects (increasing Porb) and ac- 
cretion onto the NS (decreasing Porb if Af2f > ATns) could 
have affected the orbital evolution since the SnQ However, 
we are not able to trace back the detailed evolution since 
the SN. 



^ We tried to determine Porb from the All Sky Monitor XTE data 
of CirX-1 over the last 3 years without any success. This is at- 
tributed to the poor timing resolution of the data near periastron. 
(However Fender 1997, found some evidence that the quadratic 
ephemeris is better than the linear, implying Porb 7^ 0). 
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Figure 1. The distribution of various pre- and post-SN binary parameters obtained by Monte Carlo simulations in order to construct 
systems similar to CirX-1. We assumed Aff,jg = 1.3Mq, M2 > I.OMq and 2.5 < Myi^/Mq < 6.5. We chose flat trial distributions of 
Mho) Af2i ^0 and w, and assumed an isotropic distribution of kick directions. The flgure shows the characteristics of 10* binaries (out of 
~3.4 X 10* trial systems) which passed our selection criteria for resembling CirX-1 — see text for discussions. 
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Table 2. The mean and minimum kick velocity imparted to the neutron star as a function of the assumed systemic velocity of CirX-1. 



We assume c(=5.5 kpc and 



=25' for G321. 9-0.3. v^^^ 



^rad 



and v± = <i(50/isNR- All velocities are in km s - see also Fig. 2. 
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fmin 






430 


740 


500 


no association with G321.9-0.3 and s;^ = is assumed (devils advocate's 


solution for lowest possible kick) 


430 


740 


520 


no association with G321.9-0.3, no shell impact effects included and vj_ - 


= 


580 


840 


760 


association with G321.9-0.3 and tsNR=10^ yr 




580 


850 


680 


association with G321.9-0.3, tsNR=10^ yr and -Ecjcct = 5.0 X 10^^ erg. 




700 


940 


920 


association with G321.9-0.3 and tsNR=70 000 yr 





The results of our simulations are shown in Fig. 1. We now 
summarize what can be learned from these calculations. 

The top panel shows Usys as a function of Mn, and 
^orb as a function of eccentricity. There is 76 % probability 
that 430 < Ways < 500 km s~^ and only 4.5% probability 
that «sys > 580 km s~^ if nothing is assumed to be known 
about the transverse velocity. The most probable present 
mass of the companion star, M^i strongly decreases with 
increasing values of M2i- If M2 > 1.0 M0 we find there is 
only 30% and 5.6% probability that hhi > 2.0 Mq and 
3.0 M0, respectively. We therefore conclude that CirX-1 is 
most likely to be an LMXB. An absolute upper limit for M2i 
is 4.6 Mq. 

In the right-hand panel we see that we can constrain e = 

0. 94 ± 0.04. The reason of this high expected eccentricity is 
that (as we have argued) the companion must be unevolved 
and therefore has to pass by very close to the NS in order to 
transfer mass at each passage. (Also the extreme velocity of 
CirX-1 infers that the SN nearly disrupted the system and 
thus e is expected to be large). 

In the second panel we plotted in the left-hand side the 
distribution of the allowed masses of the two stellar com- 
ponents prior to the SN. This, somewhat surprising plot, 
shows that we can not constrain the combination of these 
two pre-SN masses very well. All the combinations of masses 
shown can result in a post-SN binary which resembles Cir X- 

1. This result reflects the enormous amount of freedom due 
to the unknown direction and magnitude of the kick im- 
parted to the newborn neutron star. It should be noted that 
while the probability distribution of M2 is decreasing with 
M2, the most likely mass of the exploding He-star increases 
monotonically with Mhc- It is seen that the distribution of 
M2i is slightly shifted towards lower masses compared to 
the distribution of Al2. This refiects the mass loss from the 
companion star due to the shell impact. This effect is shown 
in the right-hand side of the second panel, where we have 
plotted the fraction of material lost from the companion 
star immediately after the SN as a function of the pre-SN 
separation, ro. We see that we expect 5-12% of the mass 
of the companion star to have been lost as a result of the 
shell impact. This is significant and future observations of 
the supernova remnant might be able to give evidence for 
this scenario. It is seen that the effect of the shell impact 
decreases with ro as expected, see e.g. equation (4). 

In the third panel we plotted ro as a function of A/2, 
and Porb as a function of Mua- The lower limit to ro as a 
function of M2 is given by the requirement that the com- 
panion star must be able to fit inside its pre-SN Roche-lobe, 
cf. equation (9). It can be seen that it is difficult to produce 



a system like Cir X-1 if ro is larger than 9 Rq (if ro is larger 
it becomes very difficult to keep the system bound after the 
SN and still require a large space velocity). Such a tight 
pre-SN binary is most likely to originate from a common 
envelope evolution. This fact also yields some preference for 
a low value of M2, since common envelopes are most likely to 
form for extreme mass ratios between its stellar constituents 
(e.g. Iben & Livio 1993). 

In the bottom panel we plotted Vsys as a function of the 
kick velocity (actually speed) and the distribution of allowed 
kick angles, in the left- and right-hand side, respectively. 
Cumulative probability curves for the kick velocity, w axe 
presented in Fig. 2. An important result is that under no 
circumstances, for any pre-SN parameters, could Cir X-1 be 
formed in a symmetric SN. In recent papers Iben & Tutukov 
(1996,1998) claimed that there is no proof for asymmetric 
SN in nature. Here our simulations show that the remark- 
ably high radial velocity observed in Cir X-1 can only be 
explained by a substantial asymmetry in the SN explosion. 
For an average set of pre-SN parameters presented in Fig. 1 
the probability for a binary to survive such a large kick is 
less than 10 %. The allowed parameter space of the kick an- 
gles show that the kick must have been directed backwards 
(i9 > 110°) with respect to the pre-SN orbital motion of the 
exploding He-star. This makes sense since a backward kick 
strongly increases the probability of a binary to stay bound 
(Flannery & van den Heuvel 1975; Hills 1983). Similarly, we 
also expect the NS to be kicked in a direction close to the 
orbital plane (<^ ~ 0° or <^ ~ 180°) in order for the system 
to minimize its orbital energy and avoid disruption. How- 
ever, if the resulting velocity vector of the NS relative to the 
companion star, just after the shell decouples gravitationally 
from the system {uq = v + w — Vim; Tauris & Takens 1998) 
is directed exactly towards the companion, then the NS is 
shot directly into the envelope of the companion star and 
the binary is assumed to merge since a(l — e) < R2f- This 
explains the nice 'shadow image' of the companion star seen 
in the plot. Hence, in the case of Cir X-1 the NS must have 
been kicked in a direction close to the companion star, but 
without hitting it, in order for the binary to survive. It is in- 
teresting to notice that the radio image of Cir X-1 (Stewart 
et al. 1993) seems to indicate that the jet is perpendicular 
to the motion of the system. Since this jet is likely to be per- 
pendicular to the inner accretion disk and the orbital plane, 
the motion of the binary seems to be in a direction along its 
orbital plane. This supports the requirement that the kick 
was directed near the pre-SN orbital plane. 

In Table 2 we list the mean and minimum values of 
the kick velocity needed to produce Cir X-1 as a function 
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Figure 2. Curves representing tlie probability that the kick 
imparted to the NS in CirX-1 was less than the value given 
on the X-axis. The full, dashed and dashed-dotted line corre- 
spond to Vbjb = 430,580 and 700 km s~^, respectively. The full 
gray line assumes no shell impact on the companion star and 
Vsys = 430 km s~l; the dotted gray line includes shell impact 



with E, 



eject 



5.0 X 10 erg and assumes Vs- 



580 km s" 



of Vs- 



It is interesting to notice that the minimum kick 
velocity needed to explain the motion of CirX-1 is smaller 
when the effects of the shell impact are included. This is 
explained by the fact that the momentum from the incident 
shell ejecta blows the companion star in the same direction 
as the NS is kicked, since the majority of NS are kicked in the 
direction of the companion in order to survive (as mentioned 
above). Hence a smaller kick is needed to produce a high 
runaway velocity, and the probability of surviving the SN is 
also increased. 



3.1 Observations of the periastron passage 

We have shown for the condition of Roche-lobe overflow at 
periastron, that CirX-1 is expected to have a highly eccen- 
tric orbit (e = 0.94 ± 0.04). Therefore the companion star 
passes by the NS at periastron in a very short interval of 
time. We find the fraction of the orbital period the compan- 
ion spends between orbital phases 9i and 62 is given by: 



1 

'2^ 



' e sin 9 



e cos 9 



+ tan 



1 + e ^ 9 

tan — 

1 - e 2 



CirX-l has an orbital period of 16.6 days. Assuming e=Q.94 
we find it only takes ~3.5 hours for the companion to move 
from 9i = 270° to 92 = 90° - see Fig. 3. The companion 
spends more than 99 % of its time between 90° < 9 < 270° . 
Hence observations must be planned very carefully in or- 
der to catch the periastron passage where e.g. the effects of 
companion heating due to X-ray irradiation are most severe. 



3.2 Possible detection of ablated material from 
the shell impact and dependence on Reject 

Our calculations show that we expect ~10% of the mass of 
the companion star to have been evaporated. Since 1.0 < 



Figure 3. Orbital geometry of CirX-1 in a reference frame fixed 
on the neutron star. We assumed M2{ = 2.0 Mq and e = 0.94. 
The periastron passage 270° < 6 < 90° only lasts ~3.5 hours. 



M2/M0 < 5.0 we estimate that ~ 0.1 - 0.5 M© of mate- 
rial in the SN nebula originates from the H-rich envelope 
of the companion star. This quantity is equivalent to ~10- 
25 % of the amount of He-rich material which was ejected 
from the collapsing He-star progenitor of the NS. We there- 
fore encourage detailed observations of the supernova rem- 
nant G321.9-0.3 to search for evidence of the shell impact. If 
-Eeject is assumed to be 5.0 x 10^^ erg then the mean fraction 
of evaporated material increases to 15 %. The constraints on 
the other parameters derived here are only slightly affected 
by a similar increase in -Eeject - e.g. the minimum required 
kick velocity is then 470 km s~^ instead of 500 km s~^. 



3.3 On the age of the SN remnant G321.9-0.3 

Our simulations show that the maximum runaway velocity 
of Cir X-1 is: u^J" — 800 km s^^ . Hence the maximum trans- 
verse velocity is constrained to be: w™"" = 675 km s~^ (given 
Wrad ~ 430 km s~^). Since the transverse velocity is inversely 
proportional to the age, tsNR of G321.9-0.3 we therefore de- 
rive a minimum age of this remnant of isNR — 60 000 yr, if 
Cir X-1 is associated with it. 
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